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Abstract. We introduce a natural definition for sums of the form 

X 



when the number of terms x is a rather arbitrary real or even complex number. The 
resulting theory includes the known interpolation of the factorial by the Y function or 
Euler's little-known formula *Yln,=i h = ~ 2 In 2 . 

Many classical identities like the geometric series and the binomial theorem nicely 
extend to this more general setting. Sums with a fractional number of terms are closely 
related to special functions, in particular the Riemann and Hurwitz £ functions. A 
number of results about fractional sums can be interpreted as classical infinite sums 
or products or as limits, including identities like 

fe(-i)' 



lim 

n — >oo 



e?( 4 ™+i)„-!-»(™+D(2^)-f J] T 1 + ~ 



fc=i 

some of which seem to be new. 
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1. Introduction 

Sums of the form Ylu=i fW) are defined classically only when the number of terms x 
is a positive integer or oo. There have been certain attempts to interpolate summations 
to non-integer numbers of terms; the most famous one is probably the interpolation of 
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the factorial, which after taking logarithms can be written as 

X 

J^lnz/ = \nT(x+ 1) 



using the well-known T function which was introduced for this very purpose (compare 
Example 4). 

We propose a systematic way to extend summations to non-integer numbers of 
terms: there is a natural and essentially unique way to do this starting from the con- 
tinued summation property Yla + Yll+i = Yll- The wor ks of Euler [E, p. 97] and 
Ramanujan [B, Chapter 6, entry 4(i)] show that both of them had also looked into a 
similar direction: for example, Euler [E, p. 101] has the formula YlZ=i v ~ l = ~ 2 In 2 
(see Corollary 14). However, we are not aware of any attempts of a treatment beyond 
sporadic examples. 

In this paper, we define "fractional sums", explore the consequences of this definition 
and derive several of its properties. While we give a precise definition in Section 2, a 
motivating special case follows from continued summation with summation boundary 
at infinity: the identity 



£+ E =E 



u=l u=x+l u=l 

certainly holds for x G N. The sum on the right makes sense classically, and so does 
the middle one: Y^L x +i f( u ) = Y^L\ f( h ' + x )- We can thus use this equation to define 
the left sum for x G C. Our Definition 3 extends this simple idea to a larger class of 
functions. The general idea is to shove the ill-defined terms to +oo where they can be 
estimated precisely in the limit. In Section 5, we explore the related definition when 
the limit — oo rather than +oo is used. 

We generalize some well-known algebraic identities from an integer number of terms 
to a fractional (even complex) number of terms (Sections 3, 4 and 5). For example, the 
finite geometric series (for i6N) 

nicely generalizes to real and even complex values of x (Theorem 11). Similarly, there 
is the well-known formula for squaring a finite sum: (a + b) 2 = a 2 + b 2 + 2ab, or more 
generally for iGN: 

/ X \ 2 X 

(1) (5z a ^) =xx + 2 5z a ^ ; 

\i/=l / v=\ n<v 

in Corollary 8, we prove an analogous formula for arbitrary x E C 

While these identities express properties of our "fractional sums" , they have special 
cases which can be rewritten in terms of classical infinite sums. This way, we obtain 
several identities involving infinite sums: some of them are well known, while others 
seem to be new, and many others can be produced in a similar spirit. We give a few 
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examples in Section 7, including a proof of the formula 

to [ £ ^'»-S*.) W -inr (l + If"] = Wexp (A _ 3 C( _ 1} _ Zg>) 

which we have not found anywhere else. 

The occurrence of special functions like T or Riemann's ( function in this example 
is no coincidence; our definition of fractional sums leads in many cases quite naturally 
to special functions. Some of the relations to ( functions are discussed in Section 6. 

Some results in this paper have been announced in [MSI]. 

We should mention that the well-known Euler-Maclaurin formula can also be used 
to estimate sums with non-integer numbers of terms provided the number x of terms 
is large. In many cases, the difference between the values of Euler-Maclaurin and 
our method tends to zero as x — > oo through the reals. However, the methods are 
fundamentally different: for example, our method does not require the integrand to be 
differentiable or even continuous. 

This paper grew out of a project at the German youth science fair "Jugend forscht" 
by the first author when he was a high school student (unaware of the works of Euler 
and Ramanujan). Both of us are no experts in this area, and we apologize if references 
or due credit are missing. 

Acknowledgments. We have received encouragement, suggestions and support 
from a number of people, including Richard Askey, Mourad Ismail, Irwin Kra, John 
Milnor, Michael Stoll and the late Judita Cofman, as well as the audiences in Miinchen, 
Stony Brook and Bremen, and in Copenhagen on the OPSFA conference. We would 
like to thank them all. 

2. The Fundamental Definition 

We begin by a natural definition for polynomials. 

Definition 1 (Fractional Sums for Polynomials). 

For a polynomial p : C — > C, let P be the unique polynomial with Y^l=\P{ v ) = P( n ) 
for all n G N. Then we define for every x <E C 

X 

:=P(x) . 

Moreover, for arbitrary a, b G C, we define 

b b-a+l 

^Tp(u) := p {y + a - 1) = P(b) - P(a - 1) . 

v=a v=l 

In this paper, we extend this concept to a larger class of functions as follows. 

Definition 2 (Approximate Polynomial). 
Let U C C and a G N U {— oo}. A function f: U — > C will be called a (right) approxi- 
mate polynomial of degree a if the following conditions are satisfied: 

• all u G U satisfy u + 1 G U ; 

• there exists a sequence of polynomials (p n ) n m of fixed degree a such that for 
every x G U , 

\f{n + x) — p n {n + x)\ — >0 as n — > +oo . 
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This is a semi- local condition and not too restrictive; only the behavior of f(x) as 
Re(x) — > +00 matters. For example, every /: C — > C with /(#) — > as Re(x) — > +00 is 
approximately polynomial of degree —00, and the functions f(x) — \nx and f(x) = ^fx 
on IR + are approximately polynomial of degree (i.e. approximately constant). The 
class of approximate polynomials is large enough for many interesting applications. 

Now comes our general definition of fractional sums. It uses the approximating poly- 
nomials as well as their fractional sums as defined above. After the formal definition, 
we try to motivate this definition and explain it in a number of special cases. 

Definition 3 (Fractional Sum and Product). 
An approximate polynomial f:U — > C of degree a G N U {—00} will be called right 
summable if for every a,b + 1 G U , the limit 

(n+b n 
£ p n {v) + (/(" + a - 1) - + 6)) 
u=n+a v=\ 

exists. In this case, this limit will be the definition for the fractional sum of f from a 
to b; we denote it by 

b b 
f(v) or briefly / . 

v=a 

Moreover, we can define fractional products by 

b / b 

f[f(u) :=expf£ln/(i/) 

u=a \v=a 

whenever In / is right summable. 




Remarks. 

• In the limit, n G N is always taken to be an integer. 

• The value of the sum is independent of the choice of the approximating poly- 
nomials p n . 

• If b — a G N, then the limit exists and agrees with the classical value of the 
sum. 

• If / is a polynomial, then Definition 3 is consistent with Definition 1, hence 
2 / = £ / f° r arbitrary complex summation boundaries. 

• We use the notation g for our fractional sum using "right summable functions" 
(where the undefined terms are shoved to +00); similarly, in Section 5, we 
introduce the symbol £ for "left summable functions" . The symbol is used 
for polynomials where both concepts coincide trivially. 

• It may be helpful to write down the first few cases of a explicitly (for a = 1 
and b — x); the cases a > 1 involve some calculations. 
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a = — oo 



£ /(") = £(/(") -/("+*)) 

x / n 

a = 0: £/(i/) = lim [x ■ f(n) + £ (f(u) - f(u + x) 



I/=l \ J/=l 



(7 = 1 



£ /(i/) = Jirn^ • /(») + ^±1) . ( /(n + i) _ /(„)) 

n 

+ £(/(")-/(" + *)) 

We will now try to motivate our definition. We begin with the case a = — oo, so 
that all p n = 0. In this case, we calculate the sum of / from 1 to x + n in two ways 

x+n x x+n n x+n 

£/ = £/ + £/ = £/+£/. 

1 1 x+l 1 n+1 



so we obtain 

(2) £/=£/+£/-£/- 

1 1 n+l x+l 

x+n 

If /(x) tends to as x — > oo, then it makes sense to require that — > as n — > oo; 

71+ 1 

the remaining two terms on the right hand side exactly yield the definition above. This 
motivates the definition for a = — oo. Note that in the last sum, the difference between 
upper and lower boundaries is an integer, so this sum is defined classically. In this 
heuristics, sums with a non-integer number of terms are denoted by our generalized 
symbol 

The next case is that of "approximately constant" functions /: suppose that / = In. 
In this case, the function / can be approximated by constants over regions of bounded 
diameters: for every bounded domain W and every e > 0, there is an iV G N such that 
for all n > N there is a C n G C for which \f(z + n) — C n \ < e uniformly for z EW . It 

n+x 

thus makes sense to estimate 2 f( u ) ~ x f( n ) (the approximating constants change 

v=n+\ 

with n, but the quality of the approximation improves for large n). This leads to the 
case o = of our definition. Note that exponentiating this example immediately leads 
to the interpolation of the factorial by the T function; see below. 

The general case is similar: in the right hand side of (2), the first and last sums 
always lead to the same term Y^=i(f( u + -0 — f( v + x ))> while the second sum can be 
estimated well by the exact formulas for polynomials. 

Note that "approximately polynomial of degree — oo" is a more general condition 
than f(x) — > as x — > oo; similarly, "approximately constant" is more general than the 
uniform condition given above for / = In: the "approximately polynomial" condition 
evaluates / only at N and at x + N, which are the only values used in our definition. 

We can now show that our definition fits in with the known interpolation of the 
factorial function. 
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Example 4 (The Extended Factorial). 
For every x G C \ {— 1, — 2, — 3 . . . }, the factorial has the following product formula 



n 

V 



Hv = lim \n x TT = T(x + 1) 



Proof. As noted above, In is approximately constant (a — 0), so our definition reads 

x / n 

(3) ^^lnz/ = lim I x In n + (ha. v — ln(z/ + x)^j 



up to an additive term in 2niZ. For every finite n, the sum is well-defined provided 
x £ C \ { — 1, —2, —3 . . . }. For v large, ln(z/) — haiv + x) = haivj (z/ + x)) is evaluated as 
the principal branch, and finitely many choices of the branch for small values of v are 
irrelevant for convergence. For Inn, we use the principal branch. 

It is readily verified that the limit in (3) exists, so the sum is well defined (up to 
finitely many additive summands 2ni, which are canceled by the subsequent exponen- 
tiation). By definition of the product in Definition 3, we get 

x / x \ / n 

v = exp I In v J = lim exp I x In n + ( In v — ha{y + x)^ 

v 



= lim ( n x "Tt 



V + X 



It is well known [C, Gauss' Formula VII. 7. 6] that this equals Y{x + 1). rj 

The following basic properties follow immediately from the definition: 

Theorem 5 (Basic Properties of Fractional Sums) . 
Fractional sums have the following properties for arbitrary a, b,c,d G C: 

b b b 

• Linearity: cJT/(z/) + d^g{y) = (cf(v) + dg(u)j , 

v=a u=a u=a 

b c c 



Continued Summation: f(y) + /(z/) = f(y) , 

b b+c 

Index Shifting: f(u + c) = ^P /(z/) , 



u=a v=a+c 

whenever two of the three fractional sums (in the last case: one of the two sums) exist. 

□ 

X — 1 X 

Note also that g / = and g / = f(x), and more generally for n G N 

x+n 



^f = f(x) + f(x + l) + ... + f(x + n). 
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3. Characterization of Fractional Sums 

The following result can be thought of as an analog to the fundamental theorem of 
calculus. We use the notation := {u ± 1 : u G U} and AS(x) := S(x) — S(x — 1). 

Lemma 6 (Summation Formula For Approximate Polynomials) . 

Let S : U — > C be an approximate polynomial of degree uGNU {— oo} ; such that G U 
(and thus N C U by definition). Then AS: U + ^ C is an approximate polynomial 

X 

of degree (7 — 1; moreover, for all x E U, the sum AS(u) exists in U and equals 
S(x)-S(0). 

Conversely, if f : U + — > C is an approximate polynomial of degree a which is right 

X 

summable, then the function S(x) := f(v) : U — > C approximately polynomial of 

v=l 

degree a + 1 and satisfies AS = f . 

Remark. The arithmetic of subtracting or adding degrees is like with differentiation 
and integration of polynomials; in particular, speaking informally, we have — 1 = — oo 
and — oo — 1 = — oo, while — oo + 1=0 (the special case — oo + 1 = — oo is possible). 

Proof. Set / := AS. Let P n be the approximating polynomials of degree a for S as 
in Definition 2, and set p n (v) := P n {v) — P n {y — 1)- For x G U + and n G N, we have 

f(n + x) = S(n + x) — S(n + x — 1) 

= P n (n + x) - P n {n + x - 1) + o(l) 
= p„(n + x) + o(l) 

as n — > oo, and / is approximately polynomial of degree a — 1. 
For all x G £7 and n G N, we obviously have 

x+n 

(4) S(n + x) = S(x) + fiy) 

v=x+l 

and hence 

x+n n 

(5) S{x)-S(0) = S{n + x)-S(n)- ^ /M + J^/M 

u=x+l v=\ 
n 

(6) 5(x)-5(0) = P„(n + x)-P n (n) + o(l) + ^(/(z/)-/(z/ + x)) . 

i/=i 

We have P„(n + x) = P„(n) + g Pn{ v ) f° r an ^ G C, so we get 

v=n+l 

n+x n 
y=n+l v=\ 

which is the definition of f(y). 



(7) S(x) - 5(0) = lim 

n— >oo 
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For the converse, we first observe that (5) holds trivially. Now we read the proof 
backwards: S is defined via (7) because / is right summable. Let p n be the approxi- 
mating polynomials of / of degree a. Then there are polynomials P n of degree a + 1 
such that 

n+x 

E Pniy) = Pn(n + x) - P n (n) , 

u=n+l 

which implies (6). Combining (5) and (6), we get 

S{n + x)- S(n) = P n (n + x) - P n (n) + o(l) . 

Starting from p n , the polynomials P n are defined only up to an additive constant, which 
we may choose so that P n (n) = S(n). The claim follows. rj 

4. Products of Fractional Sums 

In this section, we show that products like (ai + a 2 H h a n ) ■ (&i + b 2 H h b n ) 

can be multiplied out for fractional n just like for integers. 



Lemma 7 (Products of Fractional Sums). 
Let f , g : U — > C be right summable functions such that x 
an approximate polynomial. Then every x e U satisfies 



B 9iy) ) ls 



u-1 



u-1 



E/w • E^) = E /(%(") + / w • E » (*) + » w • E / (*) 



l/=l 



fc=l 



fc=l 



Proof. Let F(x) := g /(i/) and G(x) := g Then F(0) • 67(0) = and 



u=l 



A \F( x )G(x)) = F(x)G(x) - F(x - l)G(x - 1) 

= ^(x) + Fix - 1)) fax) + Gix - 1)) - F{x - l)G{x - 1) 

= f(x)g(x) + f(x)G(x - 1) + g{x)F(x - 1) . 

Since F(x)G(x) is by assumption an approximate polynomial, the first half of Lemma 6 
applies and proves the claim. rj 

Corollary 8 (Squares of Fractional Sums). 

Suppose f:U — > C is a right summable function such that x h- 
approximate polynomial. Then 



u=l 



f{u) I is an 



i/=i 



/ 2 (z,) + 2/(z,)E/W 



fc=i 



□ 



Remark. Since / is approximately polynomial, so is F := g / by Lemma 6; it is not 
automatic that F 2 is also approximately polynomial, as the example F(x) = x+sin(x)/x 
shows. 
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5. Left Summation and Binomial Series 

Classically, sums Ylv=i are defined only for integers iV with iV > 1 or iV > 0. If 
iV G N, it will be natural for us to define 

-n o 

(8) E^) ; =- E /M; 

v=l v=-N+l 

this is the only way to extend the continued summation property from Theorem 5, and 
this works for every f which is defined on {—N + 1, — N + 2, . . . , 0}. Warning: this is 

—N -N 1 

in contradiction to possible conventions like ^ = or ^ = ^ . 

v=\ v=\ v=—N 

Similarly, if x — y G — N, we set 

(9) E/M==- E /M 

u=y v=x+l 

for any / which is defined at the finitely many points {x + 1, x + 2, . . . , y — 1}. With 
this convention, Equation (8) is valid for any integer N. 

X 

Our Definitions 2 and 3 of right approximate polynomials and J? f( u ) use the 

i/=i 

behavior of / as v — > +oo; completely analogously, one can also do this for v — > — oo. 
The formal definition of left approximate polynomials is analogous to Definition 2, except 
that the limit n — > +oo is replaced by n — > — oo. 

Definition 9 (Left Fractional Sum). 
^4 left approximate polynomial f:U — > C 0/ degree a G N U {—00} m// &e called left 
summable i//or every a, 6 + 1 e (7, i/ie /zmzi 

E P» + E (/(^ + a - 1) - /(^ + &)) 

exists. In this case, this limit will be the definition for the left fractional sum of f from 
a to b; we denote it by 

b b 

or briefly ^/ . 

v=a a 

x x 

Remark. In general, the existence of the two sums £ / and ^ f is independent, and 

1 1 

if both are defined, their values can be different. 

All results about right fractional sums in Sections 3 and 4 carry over to the case of 
left fractional sums. 

Lemma 10 (Left and Right Summation). 

We have 

E^) = E^) 

v=a u=—b 

for all a, b G C for which at least one of these sums exists. 
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Proof. Define a function g(x) := f(—x). Suppose the right sum exists, so g is a 
right approximate polynomial and there are polynomials q n of fixed degree a such that 
\g(n + x) — q n (n + x)\ — > as n — > +00, for every rr. Then / is a left approximate 
polynomial with approximating polynomials q- n {—v). 

By continued summation (Theorem 5), then the definition, then changing the sign 
of the summation index, then continued summation again, and changing the sign of n, 
we obtain 

£ f(-») = - £ 9(") 

u=—b v=—a+l 

(n—b—1 n 
£ q.n(y) + E^ - °) - 9iy - b - !)) 
v=n—a+\ u=l 
-n+a— 1 —1 



= "JHM £ + £ (^( _z/ - a ) -g(-v-b- 1)) 

\i/=— n+6+1 v=-n 
/ -n-1 N 

= lim E 9n(-^) + E (d(-v-a) - g(-u - b - 1)) 



\v=— n+a u=0 

(n+b n \ b 

E + E(/(^ +«-!)- /(* + 6)) ) = E / w 

as claimed. rj 

Theorem 11 (The Geometric Series). 
For all x G C 7 we have 

q x+\ _ 1 

q u = for < q < 1, and 

o—l 
u=0 H 



~ ^ gZ + l _ 

E 9 " = q _ 1 forq>\ . 



u=0 

Proof. If q e [0, 1), we get by resolving the definition 

X CO OO x + 1 1 

E 9" = E Or 1 - O = a - E ^ = 

!,=() !/=l I/=l ^ 

as claimed. The case q > 1 is analogous. rj 

Remark. The result carries over to the case of complex q with \q\ < 1 resp. \q\ > 1, 
and the proof is the same. Care has to be taken with branch cuts even in the sum 

X 

2 (f because of the occurrence of terms q x+1 + u . It suffices to fix one branch of q x+l 

throughout the proof. 

Now we show that the Binomial series makes sense even for non-integer exponents. 
We use the general expression of binomial coefficients 

c\ c! r( c +i) 



v) ' v\{c-v)\ r(i/ + i)r(c-i/ + i) 
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Since T is meromorphic in C without zeroes and with poles exactly at non-positive 
integers, the binomial coefficient has a well-defined value in C for every c, v G C. We 
ignore the case that c is a negative integer. Then the binomial coefficient takes values 
in C, and it vanishes exactly when v or c — v are negative integers. 

Theorem 12 (The Binomial Series). 
For all c G C \ { — 1, —2, —3, • • • }, we have 

(1 + x) c = ( \x u for all x G C with \x\ < 1 and 

(1 + x ) c = ^ ( °\x v for allx G C with \x\ > 1 . 
^=o xZ/ ' 

Proof. For |x| < 1, well-known estimates imply that the summand is right summable 
with a = — oo. Therefore, we get by resolving the definition 

^>-(3 + E0'- i+ S(C)'-U e ^ 

The last binomial coefficient always vanishes and we get 

E(°y= i +E(°y=( i +*) c 

!/=0 V 7 !/=l V 7 

as claimed. For |x| > 1 we use Lemma 10 and the first part to find 

£ i !>" = e (_ c >-" = e ( c : y- = -"E r c Viy = ✓ + i 



Z/ / 14 — ' \— v — '\c—vl * — ' \V \X \ X 



and the claim follows. rj 

6. Riemann and Hurwitz Zeta Functions 

The following Dirichlet series generalizes the Riemann ( function and is known as 
the Hurwitz ( function: 



(10) {(s,x):-- 



OO 

v — 



here x is an arbitrary complex number but not a negative integer or zero. The series 
converges whenever Re(s) > I. For the definition of the powers (y + x)~ s , we use a 
branch cut at R - , and then define the complex logarithm on M~ via continuity from 
above. This way, the function x h- > C( s ) a; ) is analytic onC\ (— oo, 0]. 

For every iGC \ (— N), the Hurwitz £ function extends to a meromorphic function 
in s with a single pole at s = 1. One way to see this is via the formula 

(11) — C(s — 1, x) = — (s — 1) C(s, x) for every iGC \ (— oo, 0] 

from [SO, 64:10:1] (which is easily verified directly). 

The difference to the Riemann ( function is in the appearance of x, and in the fact 
that summation starts with v — 0; hence £(s, 1) = C( s )- 

For let [x\ denote the largest integer not exceeding x. 
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Lemma 13 (Polynomial Approximation of the Hurwitz Zeta Function). 
For every compact K C C and every s G C \ {1} ; there is a sequence of complex 
polynomials p n of degree a(s) such that 

C(s,n + z) - p n (n + z) — ► 

as n — > oo ; uniformly for z G K , with 

a(s) = 

In particular, z hh> £(s, z) is approximately polynomial of degree o~(s). 



-oo if Re(s) > 1, 
1 + LRe(-s)J i/ Re(s) < 1 and s ^ 1. 



Proof. The case Re(s) > 1 is clear with cr(s) = — oo, hence p n (z) = 0. The case that 
s G (C \ {1, 0, — 1, —2, . . . }) will be shown by induction on s — N for N e N. We may 
suppose that K is a disk with center zq. For every s ^ 1, the map x i— > Cl 5 ?^) has no 
poles for Re(x) > 0; given K, we restrict to n so all z G if satisfy Re(z + n) > 0. 

For given s, let p n be approximating polynomials for £(s, .). Let P n be polynomials 
with = — (s — l)p„ and P n (n + zq) = ((s — l,n + z ); these are approximating 
polynomials for £(s — 1, .): 

((s-l,n + z)-P n (n + z) = j (Ad/dz')C{s-l,n + z')-P' n {n + z')\ dz' 

= -(s-1) f (C(s,n + z) -p n (n + z')) dz' . 

J ZQ 

Finally, if s G {0, — 1, —2, . . . }, then £(s, z) is a polynomial in z [SO, 64:4]. rj 
Corollary 14 (Extended Riemann-Hurwitz ( Formula). 

X 

For x G C \ { — 1, —2, —3, ...} and a G C \ { — 1}, £/ie fractional sum z/ a exzsfa and 

!/=l 

satisfies 

X 

g> a = C(-a)-C(-a,*+l) . 
i/=i 

Proof. It is well known that 

(12) ((s,x+l) = ((s,x)-x- s 

for all x G C \ {-1, -2, -3, • • • } and s G C \ {1}. (For Re(s) > 1, this follows directly 
from the definition. Since both sides are holomorphic in s on C \ {1} for every x, the 
formula holds in general). 

Let S a (x) := ((-a) - ((-a, x + 1); then S a (0) = ((-a) - ((-a, 1) = and 

A5 , a (x) = S a (x) - S a (x - 1) = x a . 

Using Lemma 13, the claim follows from Lemma 6. Q 

Remark. For a = —1, we have the formula 

( 13 > E- = E(--^)=T+ S lnr(x+l)=7 + ^ + l) 
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which converges whenever x is not a negative integer. Here 7 = 0.577. . . is the Euler- 
Mascheroni constant and ip is the so-called digamma function; see [AS, 6.3.16.]. A 
special case is 



2 1 



21n2 



v 



which was noticed already by Euler [E]. 

In the following, it will be convenient to introduce the notation 



a 
X 



u=l 

By Corollary 14 and its remark, this is defined for igC\ {—1, ~ 2, —3, • • • } and a G C. 
The value x = +00 is perfectly admissible when Re(a) < — 1 and yields 

a oo = ((-a) 

with the Riemann ( function. For x G N, the function a x is known as the generalized 
harmonic series H~ a , and a x = 1 + 2 a + 3 a + . . . + x a are the first x terms of the Dirichlet 
series of £(— a) for Re(— a) > 1. 

Corollary 15 (Zeta Derivatives Identity). 
For x G C \ {-1, -2, -3, ...} ; alia G C \ {-1} and b Gff, we have 

X 

£ u a (In u) b = (-If (C (b) (-a) - C (b) (-a, x + 1)) 
where C (b) = (d/da) b (. 

Proof. This follows formally by differentiating Corollary 14 b times with respect to 
a. We omit the proof that this formal differentiation is allowed. rj 

Corollary 16 (Power Sums Up To -1/2). 
For alln G N \ {0} and a G C \ {-1}, we have 

_ i _ i 

a (-^ = £ ^ = ( 2 " C(-a) and " (^) - £ «^ = . 

Proof. This follows directly from Corollary 14 with the special values of the ( func- 
tions at the corresponding arguments. rj 

Remark. For the second identity, there is also a more direct proof which uses only 
properties of fractional sums instead those of the ( functions. By Lemma 10 and 
Continued Summation, we have for z ^ 

-1/2 -1 -1/2 -1/2 

(i4) £ ^ = £ (-"Y = - £ (-"Y = -(-ir £ v z - o 2 . 

v=\ u=l/2 v=0 v=\ 

For polynomials, left and right sum coincide trivially = £p = £p) > so ^ 

-1/2 -1/2 

z = 2n for n G N + , then § v z = — § z/ 2 , which proves the claim. 

i/=i i^=i 
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This shows that if the left and right sum both exist and are equal, there are inter- 
esting ways to manipulate those sums. Moreover, we have the following corollary: 

Corollary 17 (Power Sums and Zeros of the Zeta Function). 
Modify the definition of v z so that v z = for v — 0. Then for all z e C \ (2Z + 1) with 
2~ z 7^ 2, we have 

-1/2 -1/2 

£^ = £^ 2 <=► C(-*) = o. 

!/=l V=l 

REMARK. The definition Z := simply omits the v = 0-term in In the proof, this 
corresponds to the omission of the [y = 0)-term in (14). 

Proof. If z e C is not an odd integer, then — (— l) z ^ 1. If the two sums exist and 
are equal, then they must vanish by (14). If in addition 2~ z ^ 2, then £(— z) = by 

-1/2 -1/2 

Corollary 16. Conversely, if ((— z) = 0, then g v z = by Corollary 16 and £ v z = 

by (14). □ 

We are now going to show that series multiplication makes it possible to evaluate 
multiple fractional sums of the powers. Due to Corollary 14, the following results can 
be considered as statements about the Hurwitz ( function. 

LEMMA 18 (Double Sums of Powers). 
For arbitrary a, b G C and x e C \ {— 1, —2, —3, • • • }, we have 

X 

(15) (y -u b + b u- z/ 1 ) = ( a x) ( b x) + a+b x . 

Proof. The well-known asymptotic expansion 



Lis x) = + x — + V + , 



J ; (2j)!r( s ) 



mentioned for example in [A], shows by series multiplication that the function Z(x) := 
((—a, x + 1) • ((—b, x + 1) can be approximated by a finite linear combination of mono- 
mials (x z ) for Re(x) — > oo. Since monomials are approximately polynomial (see Corol- 
lary 14 together with Lemma 6), so must be Z. By Corollary 14, it follows that the 

function x i— > (^x^j ■ (^x^j is an approximate polynomial. 

Thus, f(x) = x a and g(x) = x b satisfy the conditions of Lemma 7. This yields 

E (£ A - ( ft) = £ v + - a E fcb + E r N 

f=l \ k=l k=l 

l/=l 

□ 



^=1 / V=l / v=\ \ k=l k=l 
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Corollary 19 (Double Power Sum). 

X 

For aeC and x eC\ {-1, -2, -3, • • • } ; we have V = a x{x + 1) - a+1 x . 

Proof. This follows from the second identity of Lemma 18 with 6 = 0, using u° — 1 
and °f = v. 

XX X 

g V + a+1 x = g (V + z/ a+1 ) = g (V • z/° + °z/ • v a ) = ( a x)(°x) + a+0 x = a x(x + 1) . 

v=\ u=l v=\ 

□ 

We can use Corollary 19 and Lemma 18 iteratively with b G N to compute arbitrary 
multiple sums of powers; a lengthy calculation gives 



x 



k=i 

x I 



jTjTr = a x( x + 1) - a+1 x 

1=1 k=l 

EEE ka = a ^Y + 3 2 x+1 )- a+lx ( x+ l) + a+2x 

rr, = -\ /—I If — 1 V ' / V / 



m=l i=l k=l 



2 



If we write ^> n+1 for the n + 1-fold iterate of the summation operator, we can write the 
general formula as 

u=0 k=l 

(this was observed by Michael Stoll). This is proved easily by induction: just show 
that by application of the difference operator A on the right hand side, n is replaced 
by n — 1, and the claim follows by Lemma 6. It is interesting to note that, after having 
defined a x for the first summation depending on a, the further iterated sums can be 
expressed in terms of a x only. 



7. Evaluation of Infinite Series, Products, and Limits 

Many infinite series can be evaluated quite intuitively with the help of fractional 
sums. In this section, we give two examples. First, consider the product 

k=l 

which has been calculated by Borwein and Dykshoorn [BD] and later by Adamchik 
[A]. We will now show that fractional sums allow to rederive a closed-form expression 
for Pi in a straightforward way. 
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Example 20 (Product by Borwein and Dykshoorn). 
For Pi : C \ {—1, —2, —3, . . .} — > C as defined in Equation (16), we have 



Pi(x) = 2~T2 



r(x + i) 



2x 



exp ( -x - 2(' ( -1, x + - ) + 2C' (-1, x + 1) - 3C'(-1) 



Proof. The function z/ i— > 2v\n (l + |) is approximately polynomial of degree a = 0, 
so by Definition 3, using lim n _ >00 2nln (l + ^) = 2x, we get 



E^(i + £) 

-i.2. + t»i»(i + f)-2(*-i)b.( 

k=i \ / \ 



lim 

n— >oo 



fc- 



Except for the summand — x, this equals the negative of the logarithm of P\ as given 
in Equation (16). Now we use the following special case of Corollary 15 

_ i 

2 

(17) £ u a In v = 2-\{-a) In 2 - (2 - 2- a )('(-a) (a G C \ {-1}) , 

together with index shifting, continued summation, Example 4 and Corollary 15: 

— - — - 

In Pi (a;) = -x -g>ln (l + -) = -x - 2 £i/ln ( ^ J 

_ i _ i 

2 2 

= — x — 2^^ v\w{v + x) + 2^^ z/lnz/ 

-x-2 £ ( z ,- a; )lnz/-^--3C / (-l) 



^=l+x 



= —x 



2 z/ In z/ + 2x In z/ 



v=\ J r x 



v=l+x 



In 2 
T2" 



3C'(-1) 



— x 



" 2 (V (-l,x + 0-C' (-1,^ + 1)) 



+2x(lnr(x + ^ - lnr(x + l^ -^-3C'(-1) 



□ 



The claim follows by exponentiation. 

Now we give an example of a limit identity that seems to be new, as far as we know. 
Example 21 (Gamma Function Product). 



lim 



2n 



k=l 



k{-l) k ' 



p V 24 2 167r 
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Remark. Using the Glaisher-Kinkelin constant A := exp — £'(— 1)), this identity 
can also be written as 

e f(4n + i) n -±-n(n+i) (27r) -f JJr 1 + - = (2e)^exp 



lim 

n^oo 



fc=l 



7C(3)\ 
16tt 2 7 



a; 



a; 



Proof. We use Lemma 7 to multiply the fractional sums lnz/ and g t>: 



w=i i/=i 

X / V 



In z/ J ■ | z/ J = — In v + ( In z/ fc J + | z/ In k 



X U=1 



k=l 



v=\ \ k=l 



Example 4 yields 



ln(x! 



x(x + 1) 



x x / 2 \ x 

-JTz/lnzy + JTlnz/fy + ^J + JTz/ln(z/!) . 

u=l v=\ ^ ' v=\ 

Specializing x = — \ and using Corollary 15 and Equation (17), we get 

_i _i _i 

1 / l\ 1 2 1 2 2 
-m7r-( — J = -^^z/ 2 lnz/ z/lnz/ + z/ln(z/!) 

V 8/ ^ l/=l ^ !/=l J/=l 



Rearranging gives 



n ln(n!) = 



to2 k(-D + k(-«)- 



n=l 



48 16 4 



8 



It is easily checked that n \— > nln(n!) is approximately polynomial of degree cr = 2. 
Resolving the definition of this fractional sum and exponentiating, while using the well- 
known identity £'(— 2) = — and Stirling's approximation of ln(n!) for large n, the 
claim follows. rj 

Many more identities can be produced in a similar spirit. Note that — \ is not the 
only interesting fractional summation boundary; other examples, resulting in interesting 
limit identities, include 



where G = Y^T=o (2fc+i) 2 denotes Catalan's constant, or 



r(i) 



explC'l-2,^ 3C(3) G 



4/ 128tt 2 4tt/ ' 



HN) ln " = exp(^ + 

n=l ^ 



7 2 7i 7r 2 In 2 2 ln 2 7r 



2 48 



+ 



where 7 = 0.577215 . . . and 71 = .072815 . . . are the Euler-Mascheroni and Stieltjes 
constants. 
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8. Questions 

The paper "On some strange summation formulas" [GIZ] contains formulas such 



as 



/,„s (— l) n /—: ^-rr 7T 2 /sin b cosft\ 

as) EH^ cosv ^ t ^ = t(^ — rj 



n=l 

and 



^ (-l) n sin ^/b 2 + vr 2 (n + 1/2) 2 _ vr sin 6 



n 71 
n=0 



The original proofs use Fourier transforms and properties of Bessel functions. It is 
tempting to prove these (and many other) identities using fractional sums. In fact, there 
is a simple argument (see [MSI, Section 9]) to derive these identities using termwise 
evaluation of fractional sums of power series: 



oo 

? 



u=a \ i=0 / i=0 \ u=a 

We leave it as an open question to establish sufficient conditions when this is justified, 
and thus to complete the proofs of (18) and (19) by fractional sums. 

A related open question is concerned with differentiation with respect to the sum- 
mation boundaries: under which conditions is the relation 

u=l v=\ 

valid with a constant c/? For example, we have (13) 

d x d x X 

— V ln(z/) = — In T(x + 1) = -7 + V - , 

where 7 ~ 0.577 ... is the Euler-Mascheroni constant. This question arises for example 
when comparing fractional sums with the Euler-Maclaurin summation formula. 



References 

[A] Victor S. Adamchik, The Multiple Gamma Function and its Application to Computation of 
Series. Ramanujan Journal 9 3 (2005), 271-288. 

[AS] Milton Abramowitz, Irene A. Stegun, Handbook of Mathematical Functions. Dover Publications 
(1974). 

[B] Bruce Berndt, Ramanujan's Notebooks, Part I. Springer Verlag, New York, 1989. 

[BD] Peter Borwein, Willem Dykshoorn, An interesting infinite product. Journal for mathematical 
Analysis and Applications 179 (1993), no. 1, 203-207. 

[C] John B. Conway: Functions of One Complex Variable I. Springer- Verlag New York (1978). 

[E] Lconhard Eulcr: Dilucidationes in capita postrema calculi mei differ cntialis de functionibus 
inexplicabilibus, 2nd ed. Commentatio 613 indicis enestroemiani, Memoires de l'academie des 
sciences de St.-Petersbourg 4 (1813), 88-119. 

[GIZ] R. William Gosper, Mourad E.H. Ismail, Ruiming Zhang: On some strange summation formu- 
las. Urinous Journal of Mathematics 37 2, (1993), 240-277. 

[M] Web page http://mathworld.wolfram.com/PiFormulas.html. 



FRACTIONAL SUMS AND EULER-LIKE IDENTITIES 



1!) 



[MSI] Markus Miiller and Dierk Schleicher: How to Add a Non-Integer Number of Terms, and How 
to Produce Unusual Infinite Summations. Journal of Computational and Applied Mathematics 
178 1-2 (2005), 347-360. 

[MS2] Markus Miiller and Dierk Schleicher: How to Add a Non-Integer Number of Terms: from Axioms 

to New Identities. Manuscript, in preparation. 
[SO] Jerome Spanier and Keith B. Oldham: An atlas of functions. Springer- Verlag Berlin, 1987. 

Institut fur Mathematik, Technische Universitat Berlin, Strasse des 17. Juni 135, 
D-10623 Berlin, Germany 

E-mail address: mueller@math.tu-berlin.de 

School of Engineering and Science, IUB: International University Bremen, Post- 
fach 750 561, D-28725 Bremen, Germany 
E-mail address: dierk@iu-bremen.de 



